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1 Introduction 

Carbon-enhanced metal-poor (CEMP) stars make up 
at least 10 % and probably as much as 20-25 % of very 
metal-poor stars with [Fe/H] < -2 l|Frebel et alj|2006l : 
iLucatello et all 120061 ). These stars have [C/Fe] > 1.0, 
and about 80 % of CEMP stars ar e also enriched i n 
s-process elements (CEMP-s stars; lAoki et al"1l2007l ). 
A likely formation scenario is pollution by mass trans- 
fer from a more massive AGB companion in a binary 
system, which has since become a white dwarf. Radial 
velocity monitoring indeed suggests that probably all 
CEMP-s stars are binaries l|Lucatello et al.ll2005aT ). 

Within the mass transfer scenario, the large pro- 
portion of CEMP-s stars requires the existence of a 
sufficient number of binary systems with primary com- 
ponents that ha ve undergone AGB nuc l eosynthesis. In 
recen t studies (jLucatello et al. 2005b; Komiva et al.l 
HmzI) it has been argued that a different initial mass 
function (IMF) is therefore needed at low metallicity, 
weighted towards intermediate-mass stars. If true, this 
in turn has important consequences for the chemical 
evolution of the halo and of other galaxies. However, 
the model calculations on which these estimates are 
based still contain many uncertainties regarding the 
evolution and nucleosynthesis of low-metallicity AGB 
stars, the efficiency of mass transfer, and the evolution 
of the surface abundances of the CEMP stars them- 
selves (e.g., see Izzard et al., this volume). 

Apart from carbon, substantial enhancements of 
nitrogen with respect to iron are common among CEMP 
stars, typically with [C/N] > 0. Detailed AGB nu- 
cleosynthesis models of low initial mass (< 2.5Mq) 
produce carbon, but do not produce nitrogen because 
it is burned during helium shell flashes. On the other 
hand, AGB models of higher mass convert the dredged- 
up carbon into nitrogen by CN cycling at the bot- 
tom of the convective envelope (hot bottom burning, 
HBB). The surface abundances of these more mas- 



sive AGB stars approach the CN equilibrium ratio 
of [C /N] sa -2. Detailed evolu tion models of AGB 
stars (Karakas & Lattanzio 2007) indicate that, while 
at solar metallicity HBB sets in at ^ 5 M© , signifi- 
cantly lower masses are needed at low metallicity (<; 
2.7 M at [Fe/H] = -2.3). One may thus expect a 
population of so-called nitrogen-enhanced metal-poor 
(NEMP) stars, with [N/Fe] > 0.5 and [C/N] < -0.5. 
Although a few examples of such stars are known, 
mostly at [Fe/H] < — 2.9, they appear to be very rare 
(| Johnson et al.ll2007n . 

In this contribution we show that the number of 
NEMP stars sets an additional constraint on possible 
changes to the IMF at low metallicity. In Sect. 2 we 
present results of our binary population synthesis sim- 
ulations, and in Sect. 3 we give our conclusions. 

2 Binary population nucleosyn- 
thesis modelling 

We have simulated populations of metal-poor halo stars 
in binary systems us ing the rapid synth etic binary nu- 
cleosy nthesis code of llzzard et al] ((20041 ) and llzzard et al.l 
(2006). The algorithm and the assumptions made in 
these simulations are briefly summarized by Izzard et 
al. (this volume). We adopt a metallicity Z — 10 -4 
([Fe/H] — —2.3) and select stars with ages between 10 
and 13.7 Gyr (roughly corresponding to the age of the 
halo) and log g < 4.0. Among this sample we designate 
as CEMP stars those with [C/Fe] > 1.0 and as NEMP 
stars those with [N/Fe] > 0.5 and [C/N] < - 0.5, fol- 
lowing the definition of I Johnson et al.l (|2007l ) . Note 
that these definitions partly overlap. We adopt a bi- 
nary fraction of 100%; for different assumptions the 
resulting numbers of CEMP and NEMP stars should 
simply be scaled with the assumed binary fraction. 

We compare our model results with t he statistics 
of the SAGA database of metal-poor stars (|Suda et al] 
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Table 1: Number fractions of CEMP and NEMP 



;ars amonj 
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2008). We selected 375 stars from the database in a 
metallicity range [Fe/H] = -2.3 ± 0.5 and log g < 4.0. 
Of these, 296 have a C abundance measurement and 69 
classify as CEMP stars, yielding a CEMP fraction of 
18-23 %, while 88 % of these classify as CEMP-s (de- 
fined as having [Ba/Fe] > 0.5). Only one star classifies 
as a NEMP star, giving a very small nominal NEMP 
fraction of ~0.3 % and a ratio of NEMP to CEMP stars 
of ~0.015. If we consider an extended metallicity range 
—4 < [Fe/H] < —2 in order to improve the number 
statistics, we find 7 NEMP stars and 78 CEMP stars, 
giving a NEMP fraction of about 1.5% and a NEMP 
to CEMP ratio of 0.09, which we regard as an upper 
limit to the NEMP/CEMP ratio at [Fe/H] w -2.3. 
We note that the ratio of NEMP to CEMP stars is 
a better constraint on our model predictions, because 
it is independent of various model uncertainties, most 
notably the binary fraction. 

The SAGA database is by no means a statistically 
complete sample. Bright giants are clearly overrepre- 
sented with respect to fainter giants and turnoff stars, 
compared to what is expected from the relative life- 
times of these phases. However, neither the observed 
sample nor our model results show a strong depen- 
dence of the CEMP fraction on evolution state (rep- 
resented e.g. by the effective gravity). Therefore the 
results we describe below are probably not greatly af- 
fected by this selection effect. 

In Table [T] we give the number fractions of CEMP 
stars and NEMP stars resulting from our models, for 
various assumptions regarding the physical ingredients 
(models 1 and 2) and the initial distributions of binary 
parameters (A-E). Model 1 is our default model, where 
we parameterize thi rd dredge-up ( 3 DUP) according to 
detailed models bv lKarakas et all l|2002h . yielding C- 
rich AGB stars for initial masses > 1.2Mq. In Model 
2 we assume much more efficient 3DUP in low-mass 
AGB stars, applying the modifications described in Iz- 
zard et al. (this volume). This results in almost all 
stars with initial masses 0.85-2.7 Mq becoming C-rich 
AGB stars and thus potential CEMP progenitors. 

The adopted binary parameter distributions are as 
follows: 

A. Default model: the initial primary masses Mi 
are distributed according to the solar neighbour- 



hood IMF as derived by iKroupa et al.1 l|l993h . 
the initial separations are drawn from a flat dis- 
tribution in log a (with a between 3 and 10 5 Rq) 
and the initial mass ratios from a flat distribu- 
tion in q = M 2 /Mi. 

B. Primary masses are distributed as in Model A, 
but mass ratios and orbital pe riods are drawn 
from t he distributions derived bv lDuquennov fc Mavod 
(|l99ll ) for the local population of G dwarfs, i.e. 

a log-normal period distribution with a broad 
peak at 170 years and a mass-ratio distribution 
with a broad peak at M2/M1 = 0.23. 

C. Primary masses are taken from a n alternative 
IMF f or the solar neighbourhood by Miller & Scalo 
(j 19791), described by a log-normal distribution 
with a median mass of 0.1 Mq and a width A log M = 
0.67, while separations and mass ratios are dis- 
tributed as in Model A. 

D. As Model C, but with a modified log- normal 
IMF with a median mass of 0.79 Mq and A l og M = 
0.5, as proposed bv lLucatello et ahl l|2005bl h 

E. As Model C using a log-normal IMF with a much 
larger median mass of 10 Mq and A log M — 
0.33, as proposed bv lKomiva et al.l (|2007l ). 

The models with default binary parameter distri- 
butions (1A and 2 A) are among those presented by Iz- 
zard et al. (this volume). As discussed in that paper, 
our default model (1A) fails to account for the large 
observed CEMP frequency, while Model 2A leads to an 
increase of the CEMP fraction to almost 10 %, much 
closer to but still short of the observed value. In both 
models the small overall NEMP fraction is consistent 
with the observations, but only Model 2A, with en- 
hanced dredge-up, can reproduce the observed small 
NEMP/CEMP ratio. 

Distributions B and C give some insight into the 
dependence of our results on uncertainties in the local 
(solar-neighbourhood) binary parameter distributions. 
Model IB leads to an increase by a factor of 1.5 in the 
number of CEMP stars and a factor of 2 in the num- 
ber of NEMP stars compared to M odel 1A, as the peak 
in the lDuauennov fc Mavorl l|l99ll ) period distribution 
coincides with the period range in which mass trans- 
fer is effective. On the other hand, in Model 2B the 
number of CEMP stars is somewhat smaller than in 
Model 2A because most CEMP stars now come from 
binaries with primary masses 0.85-1.2 Mq and q > 0.7, 
which is disfav oured by the mass-ra tio distribution. 
Adoption of the lMiller fc Scald (|1979| ) IMF (Model 1C 
and 2C) also gives somewhat higher CEM P and NEMP 
fractions than the IKroupa et al] (|l993t ) IMF, but the 
effect is modest. 

The results for distributions D and E show the ef- 
fect of a modified IMF. Mo del ID assumes the IMF 
that iLucatello et~all |2005b) suggest is required to re- 
produce the large CEMP fraction. The number of 
CEMP stars increases by a factor 2, but still falls short 
of the observed value. The discrepancy between our 
and Lucatello's results arises mainly because in our de- 
fault model the initial primary mass and period range 
contributing to CEMP stars are smaller than they as- 
sumed. With enhanced third dredge-up (Model 2D) 



www. publish . csiro. au /journals /pasa 



3 



the increase is only a modest 30%. Models ID and 
2D also show an increased NEMP fraction, by a fac- 
tor 4 over Model A, which is the result of a larger 
weight of intermediate-mass stars in this IMF. The 
NEMP/CEMP ratio is too large to be compatible with 
observations. This effect is muc h more extreme when 
we assume the IMF suggested bv lKomiva et all ([20071 ) 
which has a median mass of 10 Mq. Although this 
gives rise to a substantial CEMP fraction, close to be- 
ing compatible with observations without the need to 
assume enhanced dredge-up (Model IE), the CEMP 
stars are outnumbered by NEMP stars by a factor of 
two. This is clearly incompatible with the observed 
limits on the number fraction of NEMP stars. 

3 Conclusions 

Our binary population synthesis models show that the 
paucity of NEMP stars among metal-poor halo stars 
is incompatible with a strongly modified IMF at low 
metallicity, heavily weighted towards intermediate-mass 
stars. The observed limit on the ratio of NEMP to 
CEMP stars constrains allowed changes in the IMF 
to a median mass of at most ~0.8Mq, the effect of 
which on the frequency of CEMP stars is modest. A 
more promising explanation for both the ubiquity of 
CEMP-s stars and the near-absence of NEMP stars 
is that low-mass, low-metallicity AGB stars undergo 
much more efficient dredge-up of carbon than shown 
by detailed evolution models available to date (also see 
Izzard et al., this volume), perhaps as a result of pro- 
ton ingestion during the first thermal pulse (Cristallo 
et al., this volume). 
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